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Abstract.  It  has  been  conjectured,  for  any  discrete  density  function 
{p^}  on  the  integers,  that  there  exists  an  n^  such  that  the  n-fold  con- 
volution  fPj)  is  unimodal  for  all  n  >  n^  .  A  similar  conjecture  has 
been  stated  for  continuous  densities.  We  present  several  counterexamples 
to  both  of  these  conjectures. 

As  a  positive  result,  it  is  shown  for  a  discrete  density  with  a  conn  cted 
3-point  integer  support  that  its  n-fold  convolution  is  fully  unimodal  for  all 
sufficiently  large  n. 


1. 


1 .  Introduction 

The  limiting  distributions  of  properly  normalized  sums  of  independent  random 
variables  are  called  class  L  distributions,  (Gnedenko  and  Kolmogorov  19S4).  If 
the  random  variables  are  identically  distributed  then  the  limit  law  is  stable. 

The  problem  of  unimodality  of  the  class  L  distributions,  and  the  stable  laws  in 
particular,  has  only  recently  been  decided  after  numerous  false  proofs.  Yamazato 
(1978)  established  that  all  class  L  distributions  are  unimodal. 

In  view  of  the  apparent  unimodality  of  such  limiting  distributions,  A.  Renyi 
conjectured  for  a  discrete  distribution  f Pj }  on  the  integers  that  there  is  a  num¬ 
ber  n^  such  that  the  n-fold  convolution  (p^)  n  is  unimodal  for  all  n  >  n^  .  In 
a  similar  vein,  P.  Medgyessy  conjectured  that  for  any  continuous  probability  density 
f  there  exists  a  number  n^  such  that  the  n-fold  convolution  f*n  is  unimodal  for 
n  >  n^  .  See  Medgyessy  (1977)  for  both  these  conjectures. 

If  these  conjectures  were  true  then  it  would  be  easy  to  deduce  the  unimodality 
of  the  stable  laws.  Both  are  in  fact  false.  We  present  two  counterexamples  to  the 
conjecture  of  Medgyessy.  The  first  is  a  bounded  and  infinitely  differentiable  den¬ 
sity  with  the  property  that  none  of  its  convolutions  is  unimodal.  We  also  exhibit 
a  density  on  [0, 2x]  which  is  continuous  and  such  that  each  of  its  convolutions  is 
nowhere  differentiable.  Since  unimodality  implies  differentiability  almost  every¬ 
where,  this  example  disproves  the  conjecture  of  Medgyessy. 

As  to  discrete  distributions,  we  construct  one  which  has  the  nonnegative  inte¬ 
gers  as  its  support  and  such  that  none  of  its  convolutions  is  unimodal. 

As  a  positive  result,  we  show  in  the  final  section  that  for  any  measure  4 
with  3-point  support  (0,1,2)  the  n-fold  convolution  4°  is  unimodal  over  the 
full  range  (0, 1, . . . , 2n) ,  for  all  sufficiently  large  n.  We  conjecture  the  corres¬ 
ponding  result  for  any  measure  with  a  finite  and  connected  support  (0,1,. ..,N). 
Already  for  N«3  this  is  an  open  problem. 
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2.  Discrete  Unimodalitv 

In  the  present  section,  we  will  be  concerned  with  probability  measure  p  on 
the  integers  TL.  Let  Pjsli({j})  denote  its  mass  at  j.  Such  a  discrete  distri¬ 
bution  is  called  unimodal  if  the  sequence  [p.  , -p.}*°°  has  exactly  one  change  of 

3+1  3-« 

sign.  A  discrete  distribution  {p^ )  is  said  to  be  strongly  unimodal  if  {p^ } *  C<Ij ) 

is  unimodal  for  any  unimodal  discrete  distribution  with  connected  lattice 

support.  As  was  shown  by  Keilson  and  Gerber  (1971),  this  happens  if  and  only  if 
2 

Pj  >  Pj  iPj+i  for  all  j  e  TL.  For  instance,  the  Poisson,  geometric  and  binomial 
densities  are  all  strongly  unimodal. 


Counterexample  I. 

Consider  integers 

(2.1)  a^  =  0  <  a^  <  Sg  <  . . .  with  ah+1/ah  *  "  * 

Let  X  be  a  random  variable  with  range  A-fa^,  a^  ,  ,  ...)  and  put  Pji  =  P[X  =  ah], 

thus,  p^  >  0  and  2^p^=l.  Let  X^  ,  X2  ,  > . .  be  independent  copies  of  X  and 

S  =X,  +X-+...+X  .  Observe  that  the  range  of  S  includes  a  *  0  as  well  as 
n  1  4  n  n  u 

arbitrarily  large  positive  integers.  Hence,  Sn  can  have  a  discrete  unimodal  dis¬ 
tribution  only  if  P[Sn  =  x]  >  0  for  all  integers  x  >  0. 

Claim  0.  Let  h  and  x  be  fixed  positive  integers  such  that  >  na^  . 

Then 


(2.2)  P[Sn  “  X1  *  0  whenever  na^  <  x  <  a^+^  . 
Moreover, 

(2.3)  P[Sn  -  nah]  -  p"  and  P[Sn  -  a^]  -  npJfV^  . 
Proof.  Suppose  P[Sn»x]  >  0.  Then  x  must  be  of  the  form 

OO  00 

x  -  E  m.a  with  m.  c  TL  and  E  m  ■  n  . 

J-0  3  3  3  J-0  3 


If  tn.  = 


=  0  for  all  1  >  h  then  x  <  a.  £  m.  *  na.  .  If  m.  >  0  for  some  j  >  h 
1  -  h  j=Q  J  “h  j 

x  >  m.a  >  a.  >  a...  .  This  proves  (2.2).  A  similar  proof  yields  (2.3). 


then  x  >  mja j  >  aj  >  a}j4.i 


To  complete  the  counterexample,  one  only  needs  to  note  that  there  exists  a 
positive  integer  x  with  na^  <  x  <  a^+^  as  soon  as  a^+^  >  na^  +  1.  For  eac** 
fixed  n  this  is  true  for  all  large  h.  Therefore,  Sq  is  not  unimodal  for  any  n. 

One  may  wonder  whether  with  n  large  might  be  unimodal  at  least  over  its 

support.  The  answer  is  negative  in  general  since  one  can  arrange  that  for  each 
large  but  fixed  n  one  has 

P<S„  '  Vl1  '  <\+l  >  P|S«  '  n*hl  ‘  ph  ' 

for  infinitely  many  h.  For  instance,  p^^/pQ  n^h^O^”  ^or  *ar*»e  **  as 
2 

soon  as  p^  =  c/ (h  + 1) . 

One  may  further  wonder  whether  assuming  P[X  =  j]  >  0  for  all  j  e  Z+  (connected 
lattice  support)  might  be  sufficient.  The  following  example  shows  that  also  here 
the  answer  is  negative. 

Counterexample  II. 

Let  Xj  ,  X2  ,  ...  and  ,  Yj  ,  . . •  be  independent  random  variables  such  that 

the  X.  are  i.i.d.  with  their  common  distribution  the  same  as  in  Counterexample  I, 

1  n  n 

while  the  Y.  are  Poisson  variables  with  mean  7k.  Put  S  -  E  X,  and  T  *  E  Y,  . 

1  n  n  i-1  1  n  i-1  1 

Let  further  Z  -X, +Y,  and  S'  -  Ez-S  +T  .  Observe  that  Z.  has  the  connected 
111  n  ,  ,  inn  1 

+  1  A 
lattice  support  TL  -{0,1,2,...}. 

Let  n  be  fixed  and  choose  h  so  large  that  a^+^  >  na^  .  We  have  from 
Claim  0  of  Counterexample  I  that 


P'S;-  Vl'11  SP|Sni"*hUh£eh 


where 


S’ '„,h'“*(p|T«‘J1:  J*  Vi'1-“h)  • 

Let  b  *b  .  denote  the  smallest  Integer  >  a.  -  1  -  na.  .  If  h  Is  sufficiently 
h  n,  h  -  h+l  n 

large  then  b.  exceeds  the  mode  [i&]  of  the  Poisson  variable  T  ,  hence, 

b"  n 

c,  -  e_n^  (nM  b/b.!  .  On  the  other  hand, 
n  h 

'[Si-  Vl>  i  P'Wl1P|V01  ■  • 

Therefore, 


(2.4) 


P[S^ *  *h4-l  ”  1]  <  P^Sn“ah+l^  f°r  a11  large  h  » 


as  soon  as  e.  *o(p.  )  as  h  »,  which  is  the  same  as 

n  n+l 

bh 

(s ft)  /bh'.  -  °(P^+x^  88  h  *  °°  * 

For  large  h,  bn  >  j  «h+1  where  C  a^)  and  thus  {b^}  Increases  faster  than 
exponentially.  It  follows  that  (2.4)  holds  as  soon  as  fp^)  does  not  decrease 
too  fast,  for  instance,  only  at  an  algebraic  or  exponential  rate.  In  such  a 
situation,  P[S^« j)  increases  infinitely  often.  This  clearly  rules  out  unimodality. 


3.  Continuous  Unimodality 

One  can  easily  carry  over  the  examples  of  Section  2  to  the  absolutely  con¬ 
tinuous  case.  A  continuous  density  function  f  is  called  unimodal  if  there  exists 
a  value  x^  such  that  f  is  non-decreasing  over  (-»,Xq)  and  non- increasing 
over  (Xq  ,  <») .  The  following  is  an  analogue  of  the  above  Counterexamples  I  and  II. 

Counterexample  III. 

Let  (ph)“  satisfy  ph  >  0,  Zph“l  *nd  let  (e^Q  be  a8  in  Let 

f(x)  denote  the  density  function  which  is  obtained  by  distributing  the  mass  p^ 


5. 


uniformly  over  (a^  - 1/2,  a^  +  1/2).  Thus,  f(x) with  'l'  as  the 

characteristic  function  of  the  interval  (-1/2, 1/2). 
n 

Let  Sn  =  Zx£  where  X^  ,  "X.^  ,  ...  are  i.i.d.  with  common  density  f.  The 
*n  1=1  1 

density  f  of  Sn  satisfies 

(3.1)  f*n(x)  “  0  for  na^  +  n/2  <  x  <  a^+^  -  n/2  , 

(compare  the  proof  of  (2.2)).  Since  a^+^  >  na^  +  n  f°r  flll  large  h,  it  follows 
that  f^  can  never  be  unimodal. 

If  a  strictly  positive  density  is  desired  such  that  f  is  never  unimodal, 
one  may  start  with  Counterexample  II  and  afterwards  spread  the  mass  qj  =  P[Z  =  j] 
uniformly  over  the  interval  [j-1/2,  j  +  1/2]. 

Or  one  can  start  with  {X^,}  as  the  i.i.d.  sequence  of  Counterexample  I.  Let 
further  ZJL  =  Xi+Yi  with  the  X±  and  independent,  each  Y^  having  density  g. 

Thus,  the  Z±  have  density 

00 

fOO  =  £  Ph8(x-a.)  • 
h=0  n 

This  density  f  is  infinitely  differentiable  as  soon  as  each  derivative  of  g 

exists  and  is  bounded.  Let  us  take  g  as  the  standard  normal  density 
2  n 

«p(x)  =  (2n)"1/2e"x  12 .  Then  S ’  -  Ez  has  density 

i=l 

(3.2)  fn(x)  =  n”1^2  E  P[Sn*  j]v((x-j)/n1/2)  , 
n 

where  S  =  E  X,  . 
n  i«l  1 

Let  n  be  fixed  and  h  so  large  that  >  0,  where  ■  (a^+^  -na^)/2. 

It  follows  from  (2.2)  and  (3.2)  that 

VVi ' V  £  "'‘'W*1/2>  • 

On  the  other  hand,  using  (2.3), 


Wl>  £”'1/2HSn-S+1l9(0)  -  (2tn)-l/2»p;-lPh+1  • 

1/2 

One  has  <p(A^/n  ^  =  °^h+l^  88  ^  *  °°*  provided  {p^}  decreases  only  at  an 

exponential  or  algebraic  rate.  In  this  case,  ^(a^^-A^)  <  ^n^ah+l^  f°r 

large  h,  showing  that  f  is  not  unimodal. 

n 

Counterexample  IV. 

Observe  that  a  unimodal  density  function  is  necessarily  differentiable  almost 
everywhere.  In  this  example  we  exhibit  a  continuous  density  function  f  supported 
on  [0, 2jt)  such  that  the  n-fold  convolution  ffl  of  f  is  nowhere  differentiable 
for  every  n.  Actually,  Bogdanowicz  (1965)  already  showed  that  nearly  each  contin¬ 
uous  function  f  on  [0,2it]  has  this  property.  More  precisely,  in  the  space 
C[0, 2«]  with  supremum  norm  the  collection  of  f  e  C[0, 2jt]  with  each  convolution 
fQ  nowhere  differentiable  is  the  complement  of  a  set  of  first  category. 


Let  us  now  exhibit  an  explicit  density  with  this  property.  It  is  based  on  the 
following  result  due  to  Freud  (1962) .  A  short  proof  may  be  found  in  Kahane  (1964) . 

00 

Theorem  (Freud) .  Let  g(x) =  E  c.  cos(b^x),  where  {b^}  is  a  sequence  of 

k=l 

positive  integers  satisfying  Hadamard's  lacunarity  condition  b^+^/b^  >  q  >  1. 

Then  g  being  differentiable  at  some  point  Implies  that  c^  =  o(b^^)  as  k  >oo. 


To  construct  our  example,  let  f  be  a  probability  density  on  [0, 2ir]  and 

let  f  be  its  n-fold  convolution.  Note  that  f  is  carried  by  the  interval 
n  n 

[0,2irn].  Consider  further  the  essentially  finite  sum 


gn(x)  =  Z  fn(x  +  2rth)  . 
h=  - 

Since  g^(x  +  2it)  =  g^(x) ,  one  may  regard  g^  also  as  a  function  on  the  circle 
group  T  of  the  reals  modulo  2jt.  Relative  to  the  additive  group  T,  the  function 
g^  is  a  probability  density  equal  to  the  n-fold  convolution  of  g^  .  The  easiest 
way  to  calculate  g^  is  to  calculate  its  characteristic  function.  For  each  integer 


J  eimxgn(x)dx  =  {J  eimxg1(x)dx]n  =  {J  eimXf(x)dx}n  . 
T  T  R 


In  other  words. 


(3.3) 


gn00  ~  z  (Ym>  V1”* 

m 


as  soon  as 


(3.4) 


f  (x)  ~  Z  Yme”  for  0  <  x  <  2jt 

m 


Let  us  take 


f(x)  =  c.  -  Z  c,  cos  b  x  for  0  <  x  <  2n  , 
0  k=l  k 


while  f(x)  =  0,  otherwise.  Here,  the  b^  are  as  in  the  above  Theorem.  Further 

00 

Cq=  1/(2k)  and  >  0  (k  >  1)  such  that  Z  c^  =  c^  .  Thus  f  is  a  continuous 

k=l 

probability  density  with  g(0)  =  g (2ir)  =  0.  It  is  of  the  form  (3.4)  with  y o=cq» 

y.  =  y  = -c, /2  when  k  >  0,  while  y  =  0,  otherwise.  Thus  (3.3)  yields  that 
bk  "bk  ^  m 


(3.5) 


gn(x)  =  (cQ)n  +  2"n  Z  (-c^)"  cos  bkx  . 

k=l 


Suppose  f  were  unimodal.  Then  f^  is  differentiable  almost  everywhere. 
The  function  gfl(x)  restricted  to  (0, 2 ir)  is  the  superposition  of  the  finitely 
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many  translates  fn(x  +  2rrk),  (k  =  0, 1,  . . . , n-1),  and  thus  would  also  have  a  derivative 
almost  everywhere.  By  (3.5)  and  Freud's  theorem,  this  would  imply  that  (e^)n = 
o(b^)  as  k  oo. 


k  2 

Consequently,  if  we  choose  =  2  and  =  c/k  (k  >  1)  we  have  an  example 

of  a  continuous  density  f  such  that  none  of  its  convolutions  fn  is  unimodal.  (In 

fact,  no  convolution  is  anywhere  differentiable.) 

It  is  interesting  to  note  that,  for  any  e  >  0,  one  can  find  a  density  f^ 

uniformly  closer  than  e  to  the  standard  normal  density,  and  such  that  the  n-fold 
*n 

convolution  f  of  f  is  not  unimodal  and  nowhere  differentiable  for  any  n.  To 

e  t 

construct  such  an  f^  ,  we  start  with  a  random  variable  X  with  a  density  f  as 

outlined  in  Counterexample  IV.  Since  X  has  a  compact  support,  it  follows  from 

well-known  local  limit  theorems  (e.g.,  Petrov  (1975),  Theorems  7  or  15  of  Ch.  VII) 

that  the  density  h^  of  2  (X^  -  Nfi)/(<7  >/n)  is  uniformly  within  e  of  the  standard 

i=l 

normal  density  as  soon  as  N  is  sufficiently  large.  Since  N  ls  linearly  related 

to  the  previous  f^  ,  this  density  h^  is  not  unimodal  or  differentiable  and  neither 

*n 

is  any  of  its  convolutions  h^  .  Thus,  the  reasoning  behind  the  original  conjectures 
of  Renyi  and  Medgyessy  is  faulty.  The  central  limit  effect  is  much  too  weak  for  the 
property  of  exact  unimodality  of  high  convolutions. 


4.  Positive  Results  and  Conjectures 

Let  us  now  investigate  what  positive  results  can  be  obtained  concerning  the 
eventual  unimodality  of  sums  of  independent  random  variables.  Also  in  view  of  the 
counterexamples,  we  shall  restrict  our  attention  to  an  integer  valued  random  variable 
X  having  a  finite  support  A  with  Ac  {0, 1,  ...,N).  Let  Pj  =  P[X=  j]  thus 
A=  (j  e  TL-.  Pj  >  0] .  We  will  assume  that  p^  >  0  and  p^  >  0  and  further  that  the 
members  of  A  have  their  greatest  common  divisor  equal  to  1. 


9. 


Let 


n 


f(J)  =  Pt  z  X  =  j] 

n  i=i  1 

be  the  n-fold  convolution  of  { p ^ } .  The  support  of  f^  is  precisely  the  n-fold 

sum  An  =  A+A+.  ..  +  A  of  the  support  A  of  {p^}.  ®ne  bas  A^  c  {0,l,.,.,nN} 

and  0  e  A  ;  nN  e  A  .  In  order  that  f  .  be  unimodal  it  is  at  least  necessary 
n  n  n  J 

that  A  be  connected,  that  is,  A  =  {0, 1, . . .,nN} .  One  has  A  c  A  ,,  while 
n  n  ’  n  n+1 

G=U  A  is  precisely  the  semigroup  generated  by  A.  Hence,  in  order  that  f 
n  n 

be  unimodal  for  n  sufficiently  large  it  is  necessary  that  G  be  connected,  that 

is,  G  =  ZS+ =  {0, 1, 2,  . . . } .  This  rules  out  a  situation  like  A =  {0, 6, 10, 15}  since 

in  this  case  G  has  the  holes  (1,2, 3, 4, 5},  {7,8,9},  {11},  {13,14},  { 17 } ,  {19}, 

{23}  and  {29}. 

Since  G  contains  all  sufficiently  large  integers,  one  easily  shows  that, 
for  large  n,  the  support  Afl  of  f  has  only  a  few  holes  all  located  at  the  very 
beginning  and  very  end  of  Ar  .  One  might  therefore  conjecture  that,  for  n  large, 
f  (.)  is  unimodal  at  least  over  the  "solid"  part  of  A  .  However,  even  this  can 
be  shown  to  be  false.  Because  of  that,  we  will  restrict  our  attention  to  the  case 
of  a  connected  lattice  support 


A  =  {0,1,  . . . ,N}  . 


That  is,  Pj  >  0  for  j =  0, 1, ...,N  and  p ^ =  0,  otherwise.  The  results  so  far, 
as  well  as  certain  numerical  calculations,  lead  us  to  make  the  following  conjecture. 

Conjecture.  Let  {Pq  «  »•’•»  PN?  be  a  finite  discrete  distribution  with 

>  0  (i =  0, 1, . . . ,N) .  Then  the  n-fold  convolution  f^C.)  °f  { p ^ 3  is  unimodal 
for  all  sufficiently  large  n. 


10. 


If  N=1  then  this  convolution  is  a  binomial  law  and  thus  always  unimodal. 
Nothing  seems  to  be  known  for  the  case  N  >  2,  not  even  in  the  symmetric  case. 

The  conjecture,  if  true,  would  have  further  consequences.  For  instance,  apply¬ 
ing  it  to  the  discrete  density  p^c^p^t^  with  t  >  0  fixed,  it  would  follow  that 
also  f^(j)  =  f ^ ( j ) t becomes  unimodal  in  j  for  n>ng(t). 

For  density  functions  we  make  an  analogous  conjecture:  if  f  is  an  analytic 
density  on  the  finite  interval  (a,b)  with  only  finitely  many  modes,  then  for  n 
sufficiently  large,  the  n-fold  convolution  f  of  f  is  unimodal  over  the  full 
range  (na,nb). 

The  major  result  of  the  present  section  is  that  the  first  conjecture  is  true 
when  N  =  2 . 


Theorem.  Let  X^  (i=l,2, 
with  connected  lattice  support 
sufficiently  large. 


. ..)  be  a  sequence  of  i.i.d.  discrete  random  variables 

n 

{0,1,2}.  Then  S  “  E  X,  is  unimodal  for  all  n 
"  i=l 


Proof.  Let  p^  =  P[X  =  i]  thus  p^  >  0  for  i=0,l,2  and  p^=0,  otherwise. 
Let  fn(j)  =P[Sn=  j]  thus 

(4.1)  ffl(j)  >  ®  f°r  3  =  — ,2n  ;  fn(j)  =  0  *  otherwise  . 


Moreover, 

(4.2) 


fn+1(j>  =  £  Ptf  (l-D 

n+!  i=Q  i  n 


Let  us  introduce  the  ratios 


(4.3)  Pn(j)  -  fn(j+l)/fn(j)  , 

letting  pn ( j )  for  j  <  0  and  Pn(j)  =  0  for  j  >  2n.  Note  that  unimodality 


11. 


of  S  above  an  integer  mode  m.(n)  is  equivalent  to  p  (j)  >  1  for  j  <  mrt(n) 
n  u  n  u 

together  with  Pn(j)  <  1  for  J  >  m^n). 

Lemma  1 .  We  have  Pn(j+2)  <  PR(j)  f°r  all  j- 

Proof.  We  shall  proceed  by  induction.  The  stated  result  is  equivalent  to 

(4.4)  fn(h)fn(h  +  3)  <  fn(h  +  l)fn(h.+  2)  for  all  h. 

Since  f  (h)  =  0  if  h<0  or  h>  2n,  (4.4)  is  obviously  true  if  h  <  -1  or 

n  — 

h  >  2n  -  2,  and  therefore  for  n=  1. 

In  view  of  (4.2),  inequality  (4.4)  with  n  resplaced  by  n  +  1  is  equivalent 
to  showing  for  all  j 

2  2  2  2 

(4.5)  S  P  f  (j-r)  Z  p  f  (j +3-s)  <  E  p  f  (J  +  l-r)  E  P  f  (j+2-s)  . 

r=0  r  n  s=0  S  n  r=0  r  n  s=0  s  n 


We  must  show  that  (4.4)  implies  (4.5).  Rearranging  terms,  we  see  that  (4.5)  is 
equivalent  to 


2  2 

(4.6)  0<  E  P^[fri(j+l-r)fn(j+2-r)-fn(j-r)fii(j+3-r)] 

r=0 


+  .E 


KS 


p  p  F, 
r*s  j- 


r,  j-s 


where 


(4.7)  Fj>k  =  fn(j+l)fn(k  +  2)  +fn(j  +  2)fn(k  +  l)-fn(j)fn(k+3)-fn(j+3)fn(k)  . 

Applying  (4.4)  with  h=  j-r,  we  see  that  the  first  sum  in  (4.6)  is  non-negative. 
Hence,  it  suffices  to  show  that  >  0,  F^  ^ ^  >  0  and  ^j-1  j-2  — 

From  (4.7),  condition  F ^  ^  >  0  reduces  to 

f*(j+D  >  fn(j-l)fn<j  +  3>  . 


One  may  as  well  assume  that  1  <  j  <  2n  -  3  so  that  fn(j)fn(j+2)  >  0.  Now, 
applying  (4.4)  with  h=j  and  h  -  j  -  1  and  then  multiplying  the  results,  one 


obtains 
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fn(j”1)fn(j+2)fn(j)fn(j  +  3)  ^  fn(j)fn(j  +  1)fn(j  +  1)fn(j  +  2) 


Dividing  by  £n ( j ) ( j  +2),  one  obtains  the  desired  result.  Replacing  j  by  j  -  1, 

one  also  has  F,  1  _  >  0. 

From  (4.7),  condition  F,  ,  _  >  0  reduces  to 

J  9  J  “ 


fn(J'2)fn(j+3)  ^  fn(j"1)fn(j+2)  ' 

One  may  as  well  assume  that  2  <  j  <  2n  -  3  so  that  fn(j)fn(j+l)  >  0.  Applying 
(4.4)  with  h=  j-2  and  h=  j  and  then  multiplying,  one  obtains  that 

fn(J_2)fn(j+1)fn(j)fn(j+3)  <  +1)fn<J +2>  ' 

Dividing  by  fn(j)fR(j  + 1) ,  one  obtains  the  desired  result.  This  completes  the  proof 
of  Lemma  1. 


Remark .  Lemma  1  is  related  to  the  paper  "A  Hurwitz  matrix  is  totally  positive", 
by  J.H.B.  Kemperman,  which  has  been  submitted  for  publication. 

It  follows  from  Lemma  1  that  PR(J)  is  mondtonically  decreasing  if  j  runs 
through  the  even  integers  and  also  when  j  runs  through  the  odd  integers. 

As  to  the  Theorem,  it  suffices  to  prove  that,  for  all  sufficiently  large  n, 
there  exists  an  integer  k  =  k(n)  such  that 

(4.8)  Pn(k-2)  >  1  ;  pn(k-l)  >  1  ;  p^QO  <  1  ;  pn(k  +  l)  <  1  . 

For  afterwards  we  have  from  Lemma  1  that  pQ(J)  >  1  for  all  j  <  k-1  and 
Pn(j)  <  1  for  all  j  >  k,  implying  that  fR(.)  is  unimodal  about  k. 

In  other  words,  it  only  remains  to  show  that,  for  n  sufficiently  large,  the 
restriction  of  fR(.)  to  some  5 -point  set 


L 
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{k  -  2,k-l,k,k  +  l,k  +  2} 


Is  nonzero  and  unimodal  about  the  central  value  k.  We  will  do  this  by  using  a 

2 

local  limit  theorem.  In  the  sequel,  ^  =  EX,  0  =Var  X  while  Yj  denotes  the 

2 

j-th  cumulant  of  X.  Further,  B  denotes  positive  constant  with  B  >  |y3|/(2ct  )+l. 

Lemma  2  For  n  sufficiently  large,  n  >  the  restriction  of  fn(.) 

to  the  interval  |j-nn|  <  B  is  unimodal  about  one  of  the  two  integers  neighboring 

2 

the  value  n^  -y (2a  ). 

Proof.  Let  n  be  large  but  fixed  and  let  j  satisfy  |j-nu|  <  B.  Further 

put 

A  =  1/  (cr '/ n)  ;  x  =  (j  -  n^)/  (cr  «/n)  =  (j  -  np)A  . 


By  a  local  limit  theorem  for  discrete  distributions,  (Petrov  1975,  pages  207  and 
139),  one  has 

2 

o(2jm)1^2fn(j)  =  gn(x)e”X  ^  , 

where 

|Rn(x) |  <  Cn"3/2  , 


with  C  as  a  constant  independent  of  x  or  n.  Moreover, 

gn(x)  =  l  +  [(x3-3x)Y3/(6o3)]n‘1/2 

+  [(x4-6x2  +  3)Y4/(24a4)  +  <x6-15x4  +  45x2-15)  (y3)2/ (72a6) ]n-1  . 

-1/2  -x2/2  2  -2 

Since  x  =  0(n  )  and  e  =l-x  /2+0(n  ),  it  follows  from  an  easy  calcu¬ 

lation  that 


a(2«n)1/2fn(j)  *=  1  -x2/2  -xAy3/(2 a2)  +  pn”1  +  0(n-3/2)  , 


14. 


t  2  6 

where  P  =  yj  (8c  )  -  5  (y^)  /  (24  a  )  is  a  constant. 

Replacing  j  by  j  + 1  amounts  to  replacing  x  by  x  +  A .  Therefore, 

(2,m)1/2[fn(j  f  l)-fn(j)]  -  -4(x+A/2+A  y3/(2<72)  +  0(n-1))  . 

Here,  x  =  (j-nu)A.  It  follows  that 

fnO  +  l)  <  fn(j)  as  soon  as  j  >.  n^  -  1/2  - y^/ (2<j2)  +  0(n-1^2)  . 

Similarly, 

'  2  -1/9 

f  (j  +  1)  >  f  (j)  as  soon  as  j  <  n^i  -  1/2  - y- (2a  )  -  0(n  '  )  . 

n  n  j 

We  conclude  that  the  restriction  of  fn(.)  to  the  interval  | j  — np.  |  <  B  is  strictly 

2 

positive  and  unimodal  about  the  integer  k(n)  closest  to  the  value  nu-y3/(2o  ), 
with  one  exception.  Namely,  there  is  a  constant  K  >  0  such  that  for  integers  n 
satisfying 

|nu-Y3/(2a2)  -  1/2  - k(n) |  <  Kn"1/Z  , 

with  k(n)  as  a  (unique)  integer,  one  can  only  say  that  the  above  restriction  is 
unimodal  about  one  of  the  two  values  k  (n)  or  k (n)  + 1 . 

Remark .  It  should  be  noted  that,  under  mild  side  conditions,  the  Theorem  and 
its  proof  carry  over  to  the  case  of  independent  random  variables  X^  ,  Xj  ,  • • •  which 
are  not  necessarily  identically  distributed,  each  having  either  (0,1)  or  (0,1,2) 
as  its  support.  In  fact,  all  that  is  needed  is  that  for  n  sufficiently  large 
there  exists  an  integer  k*k(n)  such  that  the  distribution  of  Sn  *  X^  +. .  .  +  XR 
restricted  to  the  5-point  set  (k-2,k-l,k,k  + l,k  + 2}  is  unimidal  about  k.  Con¬ 
ditions  for  this  may  be  derived  from  local  limit. theorems  for  sums  of  independent, 
non-identically  distributed  random  variables  (cf.  Petrov  (1975)). 
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It  has  been  conjectured,  for  any  discrete  density  function  onthe 
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unimodal  for  all  n  A  similar  conjecture  has  been  stated  for  continuous 
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